Effect of transition-metal substitution in iron-based superconductors 
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We study theoretically the effect of transition-metal (TM) substitution in iron-based supercon- 
ductors through treating all of the TM ions as randomly distributed impurities. The extra electrons 
from TM elements are localized at the impurity sites. In the mean time the chemical potential 
shifts upon substitution. The phase diagram is mapped out and it seems that the TM elements can 
act as effective dopants. The local density of states (LDOS) is calculated and the bottom becomes 
V-shaped as the impurity concentration increases. The LDOS at the Fermi energy p(uj = 0) is finite 
and reaches the minimum at the optimal doping level. Our results are in good agreement with the 
scanning tunneling microscopy experiments. 
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The iron-based superconductors have been studied in- 
tensively since their discovery^. It is widely believed that 
the main physics in this family of materials is with the 
iron-arsenic planes. Superconductivity can be realized by 
doping either holes or electrons into the system. 

The transition-metal (TM) elements (such as Co or Ni) 
are widely used to achieve superconductivity by substi- 
tuting the iron ions 2 - - — . However, this substitution is sig- 
nificantly different from other doped materials, namely, 
the TM ions enter the conducting planes and may also act 
as the impurities. The competition of the doping effect 
and impurity effect in the TM-doped materials is of great 
interest. The scattering effect induced by the impurity 
has been studied intensively. It has been proposed to ac- 
count for many unusual physical properties^—. Since it 
is expected that all of the TM-ions enter the conducting 
planes, the impurity concentration is expected to equal 
to the doping density. While so far little attention has 
been paid to this issue when studying the impurity ef- 
fect, which may account for some unusual experimental 
observations of the TM-doped compound. 

Recently, the TM substitution effect has attracted 
broad interest - — . For Co(Ni) substitutes, it was re- 
ported that the extra electrons are concentrated at the 
Co or Ni sites based on the first-principle density func- 
tional method 1 ^—. Very recently, based on the X-ray 
absorption experiment it was also indicated that the elec- 
tronic occupations for iron sites keep constant as Cobalt 
density changes in the BaFe 2 _2a;Co2 a ;As2 compound 19 . 
These results lead to the fundamental question, namely, 
whether the TM substitutes provide doping carriers to 
the system. On the other hand, the numerical calcula- 
tion based on the first principle calculation revealed that 
the Fermi energy shifts upon TM substitution although 
the extra electrons are localized around the impurity 
sitei£~— . The angle-resolved photoemission spectroscopy 
(ARPES) experiments on BaFe 2 _2a;Co2a;As2^ i 2i have 
indicated the evolution of the Fermi surface with Co- 
substitution, which also seems to propose that the Cobalt 



atoms could be treated as dopants. Thus at this stage the 
TM substitution effect is non-trivial and rich in physics. 
Studying this issue theoretically may provide insightful 
hints to clarify the mechanism of superconductivity. 

In this Letter, we study the Cobalt doped material and 
treat all of the Cobalt ions as the randomly distributed 
impurities. The order parameters and the Fermi energy 
are obtained self-consistently based on the two-orbital 
model and the Bogoliubov-de-Gennes (BdG) equations. 
The extra electrons are totally located at the impurity 
sites and the electron filling on Fe site keeps the same 
upon substitution. In the parent compound, the spin- 
density-wave (SDW) order is revealed due to the Fermi 
surface nesting. A rigid shift of the Fermi energy occurs 
upon substitution, which would break the Fermi surface 
nesting and suppress the SDW order. As a result, the 
superconducting (SC) order shows up. These results are 
qualitatively the same with the previous ones in the clean 
system 2 ^. Thus we suggest that the Cobalt substitutes 
should be treated as the effective dopants in spite that the 
extra electrons are localized. The local density of states 
(LDOS) is also studied and we propose that the disor- 
dered impurities are necessary to elucidate some striking 
features revealed by the scanning tunneling microscopy 
(STM) experiments. 

We start from a two-orbital model including the hop- 
ping elements, pairing term, on-site interactions and im- 



purity parte d 1 ' 22 , expressed by 



H = H, 
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The hopping term H t can be expressed by 



(i) 



H.cO-to^c^c^, (2) 



where i,j are the site indices and [i,v = 1,2 are the 
orbital indices, and to is the chemical potential. The 
second term is the paring term, which reads 



ij 
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The on-site interaction term H mt can be written 



Hmt = U (ni^nifuT + U' (n»n*)»W 

T^ificr} Oliver •> (4) 

where n^a is the density operator at the site i and orbital 
/i. The inter-orbital onsite interaction U' is taken to be 
U -2J H . 

iJimp is the impurity part of the Hamiltonian. For 
Cobalt substitutes, this term should include the poten- 
tial scattering and the magnetic part. Here we consider 
the potential scattering at the Cobalt substitutes, which 
should play dominant role based on the first-principle 
calculation-^— . Then this term is written as, 



E^L 



(5) 




The Hamiltonian can be diagonalized by solving the 
BdG equations self-consistently, 
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where the Hamiltonian Hi^j U<T is expressed by, 



(6) 



-knju + [U (n ifl a) + (U - 2J H )(nifia) 



-(U - 3Jff)(n i(iier ) + E Vs6i,i m ~ *o]%<W (7) 



The SC order parameter A^-,, and the local electron 
density (ni^) are obtained self-consistently, 



(8) 

M = E Ktl 2 /^") + E NUH 1 - ( 9 ) 

n n 

Here /(x) is the Fermi distribution function and Vi^ju is 
the pairing strength. 

The LDOS is expressed by 

pM = ^[KvMEn-u) + \vfe\*8(E n + U )], (10) 



where the delta function S(x) has been approximated by 
T/7r(x 2 + T 2 ) with the quasiparticle damping T — 0.01. 

The hopping constants used in this two-orbital model 
are expressed 



a = x,y, 



__ i + i - (-iyt 3 



(11) 
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FIG. 1: (Color online) (a-c) The intensity plots of the SC 
order parameter, magnetic order and particle number with 
the impurity concentration x = 0.03, respectively, (d) The 
two dimensional cut of the above parameters along y — 18. 
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In the following presented results, the hopping con- 
stants are chosen to be ii_4 = 1.0,0.5,-1.8,0.03. The 
chemical potential to is determined by the average elec- 
tron filling per site (n) with (n) = 2 + x, where x is the 
impurity concentration. The on-site Coulombic interac- 
tion is taken as U — 3.26 and Hund coupling Jh = 1.4. 
The pairing is chosen as next-nearest-neighbor (NNN) 
intraorbital pairing with the pairing strength V = 1.0, 
which will reproduce the s± paring symmetry 2 ^ - — . The 
numerical calculation is performed on a 24 x 24 lattice 
at the periodic boundary conditions. To calculate the 
LDOS, a 40 x 40 supercell technique is used. 

We plot the spatial distribution of the SC order 
parameter A., = |E a( ,^fti+«(ii trie magnetic order 
[Mi = \ J2^( n i^t - n vl)h an( l tne P ar ticle number m = 
Y^uc n i^<y with the scattering potential V s — —3.5 in Figs. 
l(a)-l(c), respectively. Fig. 1(d) displays the two dimen- 
sional cut of all of the parameters. As seen, in presence 
of impurities, the amplitude of the SC order parameter 
is not uniform. It is reduced at and around the impurity 
sites. The magnetic spin order is antiferromagnetic along 
the y direction and ferromagnetic along the x direction, 
which is consistent with the (7r, 0) / (tt, it) SDW in the ex- 
tended/reduced Brillouin zone. This is qualitatively the 
same with the previous theoretical calculatio n 22 ' 23 . At 
and around the impurity sites, the magnetic order Mi 
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FIG. 3: (Color online) The energy dependent LDOS on the 
impurity site, NN site and NNN site of the impurity, respec- 
tively, with the impurity concentration x — 0.12. The solid 
line represents the LDOS of bulk site which is far away from 
any impurity. 



FIG. 2: (Color online) (a) The chemical potential to and the 
average particle number per Fe ion iifc as a function of the 
Cobalt concentration x. (b) The average magnitudes of the 
SC order parameter A and magnetic order M as a function 
of the temperature with the impurity concentration x = 0.03. 
(c) The average amplitudes of the SC order parameter A and 
the magnetic order M as a function of x. (d) The calculated 
phase diagram. 

is suppressed. The particle number is about 3.0 at the 
impurity sites and recovers to 2.0 when away from the 
impurity sites. This result indicates that the electron 
filling for Fe sites does not change upon Cobalt substitu- 
tion. This is consistent with the numerical results based 
on the first principle calculation^ 5 . - — . 

The average electron density per Fe site np e and the 
chemical potential to as a function of the impurity con- 
centration x are plotted in Fig. 2(a). As shown, the 
electron density is almost unchanged (« 2.0) for all of 
the impurity concentrations we consider, indicating that 
the extra electrons of impurity atoms are indeed local- 
ized around the impurity sites. This is consistent with 
the first principle calculation^ 5 . - — . In the mean time, the 
chemical potential increases monotonically with increas- 
ing substitution level x. This result confirms the numeri- 
cal results based on the first principle calcuationi 6 - - — and 
is consistent with the ARPES experiments^ 2 ^. 

The average amplitudes of the SC order and magnetic 
order as a function of the temperature and impurity con- 
centration x are shown in Figs. 2(b) and 2(c), respec- 
tively. For the fixed impurity concentration x = 0.03, as 
seen in Fig. 2(b), both the magnetic order and the SC 
order decrease as the temperature increases and two tran- 
sition temperatures are revealed. At zero temperature, as 
seen in Fig. 2(c), the magnetic order decreases monotoni- 
cally with the increasing x and vanishes around x ~ 0.06. 
The SC order increases as x increases in the low substitu- 
tion region and reaches its maximum at x = 0.06. Then it 



decreases with further increasing x. The calculated phase 
diagram is plotted in Fig. 2(d). It can be seen that the 
magnetic order and SC order coexist in the low impu- 
rity concentration region. The magnetic order decreases 
abruptly at x = 0.06, corresponding to the quantum crit- 
ical point at this concentration. This is well consistent 
with the experimental results on BaFe2-2xCo2aAs 2 28 ' 29 . 
The SC order appears as the magnetic order is suppressed 
and the SC order reaches the maximum as the magnetic 
order disappears. All of the obtained results are qualita- 
tively consistent with the previous calculation by merely 
taking into account the doping effect^. 

We have demonstrated that the Fermi energy increases 
as impurity concentration x increases. The phase dia- 
gram and the order parameters as a function of x also 
indicate that the impurities are indeed act as the effec- 
tive " dopant" , in spite that the electron densities for Fe 
sites do not change. While the fundamental question 
still remains about how the TM substitution controls the 
phase diagram. Actually, for the TM-doped materials, 
the competition of the SDW and SC orders should play 
the essential role. The SDW instability originates from 
the Fermi surface nesting 22 . It reaches the maximum 
value for the parent compound due to the perfect nest- 
ing at zero doping. The SC order is induced by the spin 
fluctuatio n 25 - 26 and also relates to the Fermi surface nest- 
ing. While for the parent compound the SC order could 
not survive due to the suppression effect by the SDW or- 
der. Both the impurity effect and the shifts of the Fermi 
energy would suppress the SDW order and then the SC 
order shows up upon substitution. Actually, here the dis- 
order and the impurities play the important role in this 
system and to some extent they induce and enhance the 
SC order through suppressing the SDW order. 

We now turn to discuss the LDOS spectra. Shown in 
Fig. 3 is the LDOS spectra near an impurity with the 
impurity concentration x = 0.12. As seen, the spectra are 
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FIG. 4: (Color online) (a) The energy dependent LDOS of 
Ba(Fei_ x Co x )2As2 superconductor in the low substitution 
region. Inset: The energy dependent LDOS with doping 
x — 0.07 without impurity, (b) Similar to (a) but in the high 
substitution region. Inset: The solid circles are the LDOS at 
the Fermi energy p(oj = 0) versus the impurity concentration 
x, with the solid lines being guides to the eye. The open cir- 
cles and dashed line are p(uj = 0) versus the doping density 
without the impurity. 



suppressed at the impurity site and the nearest-neighbor 
(NN) sites. Two in-gap peaks are clearly seen on the 
NNN site to the impurity. These results are similar to the 
previous theoretical results for single-impurity effect^. 

Now let us study the LDOS spectra in the bulk which 
are far away from any impurity site. The average LDOS 
spectra of four sites are plotted in Figs. 4(a) and 4(b), 
respectively. The bulk LDOS spectrum without the im- 
purity scattering at the doping level x — 0.07 is plotted 



in the inset of Fig. 4(a) for a comparison. As seen, for 
the clean system the LDOS spectrum has a U-shaped 
bottom, consistent with previous theoretical results 22 . 
This is understandable because the SC order parameter is 
nodeless around the Fermi surface. While interestingly as 
the impurity concentration increases the bottom becomes 
V-shaped. This is consistent with the STM experiments 
on Ba(Fei_ x Co x )2As2 superconductors^— and the pre- 
vious theoretical studies based on the T- matrix method 6 . 
Another interesting feature of the spectrum is the finite 
value of the LDOS at the Fermi energy [p{uj — 0)]. We 
plot the zero energy LDOS as a function of the impu- 
rity concentration in the inset of Fig. 4(b). p(u = 0) 
decreases with increasing x and reaches the minimum at 
the optimal SC sample x m = 0.06. Then it increases with 
further increasing x. As seen in the inset of Fig. 4(b), 
p(co = 0) scales with | x — x rn \ . This is a striking feature 
and is well consistent with the STM experiment. For a 
comparison we also plot p(to = 0) as a function of dop- 
ing density without considering the impurity scattering. 
As seen, p(co = 0) is very low and almost doping inde- 
pendent. Our results propose that the impurity effect 
is essential when explaining some striking features from 
STM experiments. 

In summary, the TM substitution effect in iron-based 
superconductors is studied by solving the BdG equations 
self-consistently. We investigate the Cobalt substitutes 
and the Cobalt atoms are considered as randomly dis- 
tributed impurities with negative scattering potential. It 
is shown that the extra electrons are localized around 
the impurity sites, while the chemical potential increases 
with the increasing impurity concentration. These results 
are consistent with the recent first principle calculation 
and the ARPES experiments. The obtained phase di- 
agram is qualitatively consistent with the previous re- 
ports of the BaFe2-2xCo2a;As2 superconductors. It is 
found that the energy dependent LDOS spectrum is V- 
shaped, and the LDOS at the Fermi energy p(uj = 0) is 
finite and depends on the impurity concentration. These 
features are consistent with the STM experiments on 
BaFe2_2i E Co2 a ;As2 material. 
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